D-brane solitons and boojums in field theory and 
Bose-Einstein condensates 

Kenichi Kasamatsu 1 , Hiromitsu Takeuchi 2 , and Muneto Nitta 3 

1 Department of Physics, Kinki University, Higashi-Osaka, 577-8502, Japan 

2 Graduate School of Integrated Arts and Sciences, Hiroshima University, 
Kagamiyama 1-7-1, Higashi-Hiroshima 739-8521, Japan 

3 Department of Physics, and Research and Education Center for Natural Sciences, 
Keio University, Hiyoshi 4-1-1, Yokohama, Kanagawa 223-8521, Japan 

E-mail: keni chi@phys . kindai .ac.jp 
Abstract. 

In certain field theoretical models, composite solitons consisting of a domain wall 
and vortex lines attached to the wall have been referred to as D-brane solitons. 
We show that similar composite solitons can be realized in phase-separated two- 
component Bose-Einstein condensates. We discuss the similarities and differences 
between topological solitons in the Abelian-Higgs model and those in two-component 
Bose-Einstein condensates. Based on the formulation of gauge theory, we introduce the 
"boojum charge" to characterize the "D-brane soliton" in Bose-Einstein condensates. 



D-brane solitons and boojums in field theory and Bose-Einstein condensates 2 
1. Introduction 

String theory is the most promising candidate for producing a unified theory of the 
four fundamental forces of nature [I]. Before the so-called second revolution of string 
theory, five kinds of theories (typel/IIA/IIB/heterotic SO/heterotic E8) were known to 
be consistent string theories, which were formulated in perturbation in terms of a genus 
expansion of the world-sheet [2] . The discovery of the Dirichlet branes (D-branes) was a 
turning point for string theorists, showing that all perturbative string theories are related 
non-perturbatively. D-branes have been found to be non-perturbative solitonic states 
of string theory and are characterized as hypersurfaces on which open fundamental 
strings can terminate with the Dirichlet boundary condition. D-branes are the most 
fundamental tool for studying the non-perturbative dynamics of string theory, which 
would enable the solving of the problems such as the dimensionality of the universe and 
the generation number of quarks and leptons 

In some field theories, topological solitons consisting of a complex of domain 
walls and strings have been proposed as a prototype of the original D-brane in string 
theory, referred to as the "D-brane soliton". The D-brane soliton was first discussed 
by Gauntlett et al. based on the hyper-Kahler nonlinear sigma model (NLaM) [3J. A 
similar argument was extended to Abelian or non-Abelian gauge theory with finite gauge 
coupling jH El El [7] . All possible composite solitons of domain walls, vortices, monopoles, 
and instantons, including D-brane solitons in field theoretical models are reviewed in 
Ref. [8]. The correspondence between the D-brane in string theory and that in field 
theory was established by the following facts [9 J : (i) The branes should have a localized 
U(l) Nambu-Goldstone mode, which can be rewritten as the U(l) gauge field on its 
world volume. The effective description of the collective motion of the D-branes is then 
given by the Dirac-Born-Infeld (DBI) action in the low energy regime [TO]. The DBI 
action is a nonlinear action of the scalar field (corresponding to the transverse position of 
the D-brane) and the U(l) gauge field [TT] . (ii) Vortex lines attached to the domain wall 
are analogous to fundamental strings because their endpoints are electrically charged 
and identical to solitons known as "Blons" in the DBI action [T2] . 

Based on these findings, laboratory experiments of D-brane physics have been 
proposed that exploit the fact that some aspects of these field theoretical models are 
closely related to the theory of condensed matter phenomena, and thus there are several 
systems appearing in nature which can be said to admit analogs of D-branes. The 
purpose of the analogy argument is to find a way to observe at least some parts of the 
original complex theory in cosmology or high energy physics by simulating the analogous 
phenomena in condensed matter experiments. To date, D-branes have only been a 
theoretical hypothesis and it is not clear whether or not they can be used to describe 
nature or not. Therefore, it is important to simulate D-brane physics in a laboratory to 
raise the practical applicability of this hypothesis and to further develop D-brane physics. 
For example, analogues of "branes" have been discussed in an experiment on superfluid 
3 He [13] , where Brane-anti-brane collision was simulated using the boundary between 
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3 He-A and 3 He-B phases in a well controlled manner; the brane-anti-brane collision 
could explain the Big Bang in brane cosmology, after which lower-dimensional branes, 
e.g., cosmic strings, remain as relics. However, their exact correspondence to D-branes 
in string theory was not clarified in Ref. [13]. In other systems, e.g., ferromagnets, it is 
difficult to create and observe three-dimensional topological defects experimentally, so 
this kind of study has not yet been reported. 

Recently, we reported that an analogue of a D-brane can be realized in ultra-cold 
atomic Bose-Einstein condensates (BECs) and the experimental verification of various 
D-brane dynamics could be possible [151 US H7J- The domain wall in two-component 
BECs can mimic the fundamental properties of the D-brane, because the theoretical 
formulation in two-component BEC can be mapped to the NLcrM by introducing a 
pseudospin representation of the order parameter [TBI HHJ |20]. It can then be shown 
that the domain wall in a two-component BEC has a localized U(l) Nambu-Goldstone 
mode, which is related to the local U(l) gauge field on the wall as a necessary degree 
of freedom for the effective description with the DBI action of a D-brane pm [111 EE] . 
Thus, the resultant wall-vortex composite solitons could possess the characters (i) and 
(ii) described above and correspond to the D-branes in the sense of Ref. [3J. Here, 
the domain wall and the vortex can be identified as the D-brane and the fundamental 
string, respectively. Compared to the above systems, two-component BECs have a 
great advantage for exploring D-brane physics, because several experimental techniques 
to create, control, and observe the topological defects have been very well established 

In this paper, we discuss the correspondence of topological solitons between BEC 
systems and gauge theoretical models to clarify the similarities and differences. We study 
various topological solitons such as vortices, domain walls, and their complexes, i.e., D- 
brane solitons, in the Abelian-multi-Higgs model, which is a similar formulation to the 
two-component Gross-Pitaevskii (GP) model that governs the dynamics of BECs within 
the mean-field level. The GP model corresponds to the strong gauge coupling limit of 
the Abelian-Higgs model, where the gauge field is not an independent dynamical degree 
of freedom and is consistently determined by the configuration of the matter (Higgs) 
fields. We also focus on the topological structure of the wall- vortex connecting points in 
D-brane solitons. In a field theoretical model, these points form defects called "boojum" 
as the negative binding energy of vortices and a wall, and a half of the negative charge 
of a single monopole [6] [2TJ [22] . Boojums are known as point defects existing upon the 
surface of the ordered phase; the name was first introduced to physics by Mermin in the 
context of superfluid 3 He [23]. Boojums can exist in different physical systems, such as 
the interface separating A and B phases of superfluid 3 He [211 12S], liquid crystals |26j, 
the Langmuir monolayers at air- water interfaces [27], multi-component BECs with a 
spatially tuned interspecies interaction [281 EE], and high density quark matter [30J. In 
the present model, boojums can be found at the end points of vortices on the domain 
wall, at which the vortices change their character from singular to coreless type. 

This paper is organized as follows. In Sec. El we first review briefly topological 
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solitons in the Abelian-Higgs model. The analysis of this section can be applied directly 
to two-component BECs. In Sec. [3], we discuss the formulation of two-component BECs 
and discuss the similarities and differences by comparing the structure of topological 
solitons such as vortices and D-brane solitons of the field theoretical model in Sec. |2j 
We conclude this paper in Sec. HI 



2. D-brane solitons in field theoretical models 

As described in Sec. HJ similar structures of D-branes may occur in various field theories. 
In a field theory, both branes and strings arise as solitonic type objects. In this section, 
we study prototypes of topological solitons including vortices, domain walls, and D- 
brane solitons in gauge theoretical models. The analog of a D-brane in the NLcrM was 
first pointed out by Gauntlett et al, where they referred to it as a Q-kink-lump solution 
[3]. It turns out that the low energy effective theory of the collective coordinate for 
this solution can be described by the DBI action [10]. The vortices attached to the 
wall are identified as fundamental strings since their endpoints are electrically charged 
and identical to Blons in the DBI action [12J. Therefore, the NLcrM offers a simplified 
model for studying D-brane dynamics, that is instructive for studying full string theory. 
A similar discussion has subsequently been made for supersymmetric gauge theory [I]. 
For the strong coupling limit, the model is reduced to the NLcrM and the analytical 
solutions of the D-brane soliton can be obtained (6]. This gauge theoretical model is 
closely related with the two-component GP model, where two-component scalar fields 
(order parameters) represent the condensate wave function of the BECs. The analysis of 
this section can be directly applied to our problem in some limit. Importantly, from this 
formulation, we can define a suitable topological charge that characterizes the D-brane 
soliton from which we can understand the details of the connecting points, i.e., boojum, 
of the wall and the vortex. 



2.1. Vortex in Abelian-Higgs model 



We start from the simplest Abelian-Higgs model to understand the structure of the 
topological defects in the field theoretical model, which is useful in the following 
discussion. The Abelian-Higgs model is given by 



H 



d x 



2 -v 2 ) 2 



(1) 



^(E 2 + B 2 ) + |(V-zA)0| 2 + ^ 

where e is the gauge coupling (electric charge), A the Higgs scalar coupling and v 
the vacuum expectation value (VEV) of the scalar field 0, and E and B are the 
electronic and magnetic field, respectively. For the vector potential A, we have taken the 
transformation from the conventional notation as eA — > A for convenience. The model 
is identical to the Ginzburg-Landau free energy in the theory of type-II superconductor. 
The hamiltonian is invariant under the local U(l) transformation c/>(r) —» e ia ^(j)(r) and 
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Figure 1. Typical structure of the ANO vortex. 



A — > A + Va(r). The stationary point of this hamiltonian gives the equations 

(V-*A) 2 0=^(H 2 -^ 2 )0, (2) 

j = ie 2 [0*V0 - (V0*)0] + 2e 2 |0| 2 A = i [0*D0 - (D0)*0] , (3) 

where D = V — iA is the covariant derivative and the latter equation gives the electric 
current density. 

First, we give the solution of the vortex type, known as the Abrikosov-Nielsen- 
Olesen (ANO) vortex [2HE2]- By using = |0|e ie , Eq. ([3]) gives the vector potential 
as 

Without any current, the magnetic flux is quantized to be an integer (n = 0, 1, 2, • • •) 



d 2 xB z = (h di ■ A = (h di-V9 = 2vm, (5) 



as, 



where the vortex winding number is given by the first homotopy class: n G 7Ti[U(l)] = Z. 

The structure of the ANO vortex solution is well known based on the numerical 
solution of Eq. (T5]). The qualitative feature of the fields (p and B z is represented in 
Fig [TJ The Higgs field is reduced from the VEV at the vortex core, in which U(l) 
gauge symmetry is recovered. The magnetic flux is concentrated around the vortex 
core. There are two length scales over which the above two fields vary spatially. The 
gauge mass m v ~ y/2ev and the scalar mass m s ~ a/Aw give the penetration depth 
r v = m~ l ~ (y/2ev ) _1 and the coherence length r s = m" 1 ~ (a/Av) -1 , respectively. The 
ratio f3 = m 2 s /m 2 = A/2e 2 is the only free parameter in the model. For r v < r s (/3 < 1), 
the superconductor is classified as "type I". Then, the interaction between vortices is 
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attractive due to the exchange of the scalar fields and is consequently unstable under 
a magnetic field. For r v > r s ((3 > 1), the superconductor is classified as "type II". 
The exchange of the gauge fields leads to a repulsive interaction between vortices and 
is stable under a magnetic field. The case r v = r s (f3 = 1) is known as the "critical 
coupling", where there are no interactions between static vortices. 

For critical coupling (A = 2e 2 ), the equations for the topological defects can be 
simplified by taking the Bogomornyi-Prasad-Sommerfield (BPS) bound [33l 134] of the 
energy as 

H= Jd 2 x \D x <P + iD y <P\ 2 + ^{B z -e 2 {\<P\ 2 -v 2 ) 2 Y 

> -v 2 d 2 xB z = 27rnv 2 } (6) 



- v 2 I d' z xB z 



where we have integrated in the uniform z-direction. When the equality holds, we can 
obtain the energy minimum and thus the most stable configuration for a fixed vortex 
number n. The minimized condition yields a set of simple linear equation known as the 
BPS equation 

(D x + iD v )<f> = 0, B z = e 2 (\<f ) \ 2 -v 2 ). (7) 

Fields satisfying the BPS equations are automatically minima of the energy for a given 
n, and hence are guaranteed to be stable. 

2.2. Vortex in Abelian-two-Higgs model 

Next, we consider the Abelian-two-Higgs model [35]. The hamiltonian is given by 
the direct generalization of the Abelian-Higgs model where the complex scalar field is 
replaced by an SU(2) doublet $ T = $ 2 ) as 

f d 3 x\^(E 2 + B 2 ) + |(V - iA)$\ 2 + ^($t$ _ ^21 (g) 
j ize 4 

This model is similar to the gauged two-component Ginzburg-Laudau model used to 

describe unconventional superconductors. This hamiltonian is invariant under the local 

U(l) gauge transformation U(1) L : $ — > e' a ( r )$ and A — > A + Va(r), and the global 

SU(2) operation SU(2) G : $ — > e l7 ' <T $, where 7 = 711 with the positive constant 

7 G [0, 47r) and unit vector n, and a is the Pauli matrix. Because the isotropy group 

is U(1) L+G [36J, the order parameter space becomes [U(1) L x SU(2) G ]/U(1) L+G ~ S* 3 . 

The first homotopy group of the order parameter space is trivial 7Ti(5' 3 ) = (simply 

connected), so there are no topological vortex solutions. On the other hand, looking 

only at the gauged part of the symmetry, the local U(1) L symmetry actually breaks as 

in the case of the Abelian-Higgs model and thus 7r 1 [U(l) L ] = Z; the U(1) L gauge orbit 

away from the degenerate SU(2) G orbit is still expensive in terms of gradient energy, so 

it is possible that a vortex is stable. This vortex is called a semilocal vortex, defined in 

a non-topological manner [35J. 



H 
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After the symmetry breaking, there are two Nambu Goldstone bosons, one scalar 
of mass m s ~ yXv and a massive vector boson of mass m v as in the Abelian- 

Higgs case. In the special case for critical coupling (3 = A/2e 2 = 1, one can take the 
BPS bound in the same way as for Eq. fl6]): 

(D x + iD y )$ = 0, B z = e 2 (&<S>-v 2 ). (9) 

where the field <fi is just replaced with $. Thus, for the fixed topological sector, one 
can obtain the configuration satisfying the BPS equations, which is a local minimum 
of the energy and automatically stable, even though the vacuum manifold is simply 
connected. In this case, the vortex energy is independent of its core size, and thus the 
vortex can have an arbitrary size. For 0^1 the stability of the vortices depends on 
the dynamics and is controlled by the parameter /3 [37]. When the vortex energy is 
written in terms of its size, the energy monotonically increases as a function of the size 
for (3 < 1 (type I case). Then, the vortex tends to shrink to zero size, which means that 
only one component of $ survives while the other disappears, resulting in a stable ANO 
vortex. On the other hand, the energy becomes a monotonically decreasing function 
with respect to the size for > 1 (type II case). In that case, the vortex is unstable for 
expansion and its winding eventually disappears. 

For e — > oo, the electromagnetic energy (the first term of Eq. ([8])) vanishes, so that 
the vector potential A becomes non-dynamical variables, given by 
jtfy (y«t)« 

2 $t$ v ; 

In addition, for the limit A — > oo, the radial degree of freedom of $ is frozen and the 
minimum of the potential yields the VEV = |$!| 2 + |<l> 2 | 2 = f 2 . As seen below, the 
order parameter space can be reduced to S 3 /U(1) L ~ S 2 (sphere), and the corresponding 
model is also reduced to the 0(3) NLaM (CP 1 model). The order parameter can be 
written with the stereographic coordinate u by 

* = -r^=( 1 ) (11) 

and the corresponding unit isovector (pseudospin) s can be given as 

$t<r$ f u + u* u-u* 1 - |d 2 \ , , 

-«m7T5>TTTn* > ( 12 ) 



$t$ \1+ | U |2' 1 + | M |2' 1 + \ U \ 

where u = and oo corresponds to $ T = (v , 0) and $ T = (0, v), respectively, and, in 
terms of the pseudospin, it corresponds to the north and south pole, respectively, as 
shown in Fig. |2j The hamiltonian Eq. (jSJ) is reduced to 

ff =/ A Srw = j/ A E(v«.)' («) 

In this limit, the semilocal vortex is reduced to a coreless vortex, known as a lump 
(2D skyrmion) |38j; its spin configuration is shown in Fig. [3j The second homotopy 
group of S 2 is nontrivial as tt 2 {S 2 ) = 7r 2 [SU(2)/U(l)] ~ 7r 1 [U(l)] ~ Z, which gives a 2D 
skyrmion charge; the hedgehog configuration of the pseudospin can be mapped to the 
configuration of the 2D skyrmion through the stereographic projection. 
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South pole 



Figure 2. Stereographic projection from the sphere to the tangent plane at the north 
pole. 




Figure 3. Spin profile of a lump (2D skyrmion, coreless vortex). The spin orients 
upwards at the center and it continuously rotates from up to down as it moves outward 
radially. 



2.3. D-brane solitons in the massive Abelian-two-Higgs model 

In the following, we assume a critical coupling A = 2e 2 . This regime allows us to make 
an analytical treatment with the BPS equation as shown above. Our main target is to 
analyze D-brane solitons, which consist of vortices and domain walls. However, since 
the order parameter space of the model Eq. (jSJ) is S 3 , there cannot be a domain wall. 
To realize a discrete ground state configuration in space, the mass for the Higgs fields 
should be introduced. We extend the model Eq. (jSJ) to the massive Abelian-two-Higgs 
model, given by 



H = I (I'.r 

,2 



i_(E 2 + B 2 ) + |(V — *A)$| 2 + -^(VS) 2 



L2e 2V ' ,v ' 1 2e 2 

e i 

mass matrix 



+ i.($t$_ t; 2)2 + <|) t( SI _ M )2 $ j > ( 14 ) 

Here, we have introduced the neutral scalar field S, 2 x 2 identical matrix I and the 
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The role of the neutral scalar field E is to generate different vacua or, in other words, 
to determine the direction of the symmetry breaking from SU(2) G to U(1) G of the $ 
field; for example, see jU [8]. Actually, there are two ground states as ($) T = (v,0), 
(E) = Mi, and ($) T = (0,t>), (E) = M 2 . Thus, one can make a domain wall between 
these two vacua. For a single domain wall, the BPS bound is given as 

H = j dz{\D z <S> + (EI - M)$| 2 + ^[B Z - d z Y> - e 2 (|$| 2 - v 2 ) 2 ] 2 

+t w + d x J z } (16) 
with the topological charge density 

t w = v 2 djl (17) 
and the current 

J z = -$ f (SI- M)$ (18) 

After the integration, the contribution of the current vanishes and the total topological 
charge (tension) of the domain wall is 

dzt w = v 2 AM, (19) 

associated with the difference of the mass AM = Mi — M 2 between the two vacua. 

For e — > oo, both the vector potential A and the scalar field E are no longer 
independent dynamical fields and their form is determined by the minimization of the 
energy with respect to A and E, where A is given by Eq. ( TT0|) and E is 

„ $tM$ 



Mi 




2 + M 2 


$ 2 


2 






2 + 


$ 2 


2 



When we choose Mi = M/2 and M 2 = —M/2, E is identical to the ^-component of the 
pseudospin field: 

(2D 

Since the limit e — > oo means A — > oo here, the model also reduces to the massive 0(3) 
NLcrM with the relations Eqs. ( 11 ip and ( II 2p . The hamiltonian can be rewritten as 



M 


$i 


2 _ 


$ 2 


2 


2 


$i 


2 + 


$ 2 


2 



^(V Sq ) 2 + M 2 (1-, 2 ) 



(1 + |m| 2 )^ 

The exact solutions of the BPS equation can be found to be [3J 

Here, Zo represents the position of the flat domain wall (s z = 0) whose transverse shift 
causes the Nambu-Goldstone mode due to breaking of the translational invariance. The 
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(b) 



(AMy 



M 

V 1 



Mr 



(c) 



ANO vortex 



domain wall 




AM/e 2 v 2 



{ev)- 




oojum 



ANO vortex 



Figure 4. Typical examples of a domain wall and a wall-vortex complex (D-brane 
soliton). The profile of the domain wall for (a) the e — > oo limit and (b) finite e, where 
$1, $2, and £ are represented by solid, dashed, and dotted curves, respectively. Figure 
(c) represents the typical configuration of a D-brane soliton for finite e. Two vortices 
end on the domain wall from both sides. 



phase 6q corresponds to the azimuthal angle of the pseudospin s, causing the breaking of 
the global U(l) locally along the wall. By promoting these two variables to dynamical 
fields, we can construct an effective theory of the domain wall. The low-energy dynamics 
of a single domain wall is simply described by the DBI action [2111], where the local U(l) 
gauge fields living on the wall are created by the dual transformation of the localized 
zero mode of the phase 8q due to the breaking of the global U(l) in the original system. 
This is why this domain wall can be identified as an analog of a D-brane, as stated in 
Seem 

The structure of the domain wall profile at finite e depends on the dimensionless 
parameter e 2 v 2 /AM 2 jl]. A typical profile of the domain wall is shown in Fig. HI For the 
sigma model (strong gauge coupling) limit e 2 v 2 ^> AM 2 all fields interpolate between 
their vacuum values over a single length scale (AM) -1 [see Fig. IU(a)]. In the opposite 
weak gauge coupling limit e 2 v 2 <C AM 2 , the domain wall has a three-layer structure, 
as shown Fig. H^b). In the two outer layers, each field drops exponentially to zero 
with the width (et>) _1 . In the inner layer which has width AM/e 2 v 2 , the adjoint field 
E interpolates between its two vacuum values. 
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Next, we consider the composite soliton that consists of the vortices and the domain 
walls. Combining Eqs. and ffl6|) . the BPS bound under the fixed topological sector 
is given by 

H = J d 3 x{\D z $ + (ZI-M)$\ 2 + \D x <S> + iD y $\ 2 

+t w + t v + t m + V- j} (24) 

with the topological charge densities 

t w = v 2 d z Z, (25) 
t v = -v 2 B z , (26) 

t m = 4B-VS = iv-(SB), (27) 

e z e A 

and the current 

J = (i&D y ®, -i&D x <f>, -$ f (SI - M)$) . (28) 

An important result is the appearance of the "monopole charge" t m of Eq. ( 1271) . We 
can write B • VS as V ■ (BE) where V • B = holds. Note that the charge is defined by 
the projected magnetic field SB. This monopole should be distinguished from an usual 
hedgehog of S 2 , because it cannot exist in U(l) gauge theory. This point-like defect is 
expected to be located on the domain wall as the connecting point of a domain wall and 
a vortex, so that this charge may be referred to as a "boojum charge". The topological 
charge is given by integrating the topological charge density in the corresponding spatial 
dimension as 

T w = J dzt w = v 2 (M 1 - M 2 ), (29) 
T v = J d 2 xt v = 2nnv 2 , (30) 

d 3 xt m = J d 3 xV ■ (SB). (31) 

The typical configuration of a D-brane soliton is depicted in Fig. IH(c). Since there is 
only one component of $ on either side of the wall, the vortices should be identified as 
ANO vortices. The boojums are expected to be spread inside the domain wall on which 
the two vortices end from both sides, as shown in Fig. 11(c). 

In the strong gauge coupling limit, the analytic solution for the BPS wall-vortex 
soliton can also be derived as [6] 

$(w, z) = I H (w)e Mz . (32) 

Here, Hq(w) is the Moduli matrix given by 



w 



#oH= nC 1 K ■ (33) 
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In this notation, the k{ individual vortices are positioned at w — in the $j domain 
(i = 1,2). For M = diag.(M/2, — M/2), the solution can be written as 

$( "' Z) = Vl + e-™(»-*)|Z|* ( e-^-^Z 1 ' (34) 

Here, the holomorphic function Z = Z(w) represents the contribution of the coreless 
vortex (lump): 



Z W = ^ fef (35) 



This form of the solution is equivalent to that found by Gauntlett et al. [3J, where 
u(w,z) = e~ M ( z ~ Zo ) +l6o Z(w) by using the transformation Eq. (ITTj) . and they consider 
the case k 2 = k± — 1; Isozumi et al, on the other hand, consider the more general 
case including k\ = k 2 |6j. We can thus construct solutions in which an arbitrary 
number of vortices are connected to the domain wall. In the BPS solution, the energy 
is independent of the vortex positions Wj ' on the domain wall; in other words, there is 
no static interaction between vortices. The structure of a D-brane soliton with a simple 
wall-vortex configuration is shown in Fig. |5j Here, the $i and $2 fields are positioned 
in the z > and z < regions, respectively. In Fig. Et^a), a single vortex (lump) exists 
in z < ($2 domain), given by Z = w. 

The edge of the vortex attaches to the wall, causing it to bend logarithmically 
as z = log|w|/M. Figure GD^b) shows a solution in which both fields have one vortex 
connected to the wall, corresponding to Z = (w — x ^/(w-xjp). In this case, the wall 
becomes asymptotically flat due to the balance of the vortex tension. 

As shown above, the connecting points of the wall and vortices can give rise to 
point defects known as boojums (6j ETJ [22]. The presence of a boojum is confirmed by 
calculating the boojum charge by Eq. ( 13T|) . Because the domain wall is logarithmically 
bent when a single vortex is attached to one side, this situation does not give a clear 
integral region to calculate the topological charge. Let us consider that two vortices with 
winding number n = 1 are attached to both sides of the wall. Then, we do not need to 
take care of the integral region of Eq. (13T1) because of the asymptotically flat wall. The 
charge T m in this configuration, i.e., the charge of two boojums, can be calculated as 

T m = 4 / ^/ : VV ■ (HB) 



e 2 
1 

^2 



d 2 xB 7 



d 2 xB 7 



_^( Ml -M 2 ). (36) 



e 2 



It turns out that boojums have negative monopole charge and negative energy, despite 
the signs of the charges of the vortices and domain walls [21J. Since one monopole 




Figure 5. Typical configuration of a D-brane soliton in the field theoretical model 
under the sigma model limit. The left pictures show the isosurface of the domain wall 
u = 1 (s z = 0). (a) A single vortex is attached to the domain wall from one side, (b) 
Two vortices are attached to the wall from both sides. The spin configuration in three 
z =const planes, chosen as z > 0, z — 0, and z < 0, are also depicted on the right. 
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Figure 6. Schematic illustration of the wall- vortex soliton configuration in two- 
component BECs. The two-component BECs (z > 0) and ^2 (z < 0) are separated 
by a domain wall near the z — plane. The domain wall can be defined as a plane 
on which both components have the same amplitude. We assume that vortex lines, 
nucleated by rotation around the z-axis, are perpendicular to the wall. 

is characterized by T m = |f (Mi — M 2 ) > [39J, one boojum has a half of the single 
negative monopole charge. It has been proposed that the boojum plays a role in binding 
the wall and the vortex, because the energy of this composite soliton is smaller than the 
independent sum of the energies of the wall and the vortices. 

3. D-brane solitons in two-component BECs 

Based on the above analysis, we consider a system of cold atomic two-component BECs 
and study the detailed properties of the wall-vortex complex (D-brane soliton) in this 
system. The system under consideration is schematically shown in Fig. [6] Two- 
component BECs have been realized by using a mixture of atoms with two hyperfine 
states of 87 Rb [101 SH S2] or a mixture of two different species of atoms such as 87 Rb- 
41 K pi E], 85 Rb- 87 Rb pE], and 87 Rb- 133 Cs [38]. The experiments in Refs.pEJ SH] 
demonstrated that the miscibility and immiscibility can be controlled by tuning the 
atom-atom interaction via Feshbach resonances. Here, the domain wall is referred to 
as a boundary of phase-separated two-component BECs and is well-defined as a plane 
on which both components have the same amplitudes [HI SHI SHI EH [51]. The vortices 
can be arranged by applying a rotation to the BEC around the z-axis or imprinting an 
appropriate phase profile to the condensate with engineered atom-laser coupling [52] . 
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3.1. Gross-Pitaevskii model 

Two- component BECs are described by the order parameters 

* = i ;/ 1 = i ^ e 2 i , (37) 




p2e 

which are condensate wave functions with density pj and phase 9j (j = 1,2). They 
are typically confined in a trapping potential Vj. The energy functional of the two- 
component Gross-Pitaevskii (GP) model is given by 

w_ 

.2m a 



ff[*i,* 2 ] = / d 3 x{j2[ 

J 3=1,2 



^ .2m.- -A T 2 



h 



+ (V J - f i j )\^ j 



2 



+ ¥l*j| 4 



+ <7i2|*i| 2 |^2| 2 }. (38) 



2 

Here, rrij is the particle mass of the j-th component and fij is its chemical potential. 
The coefficients gxi, g 22 , and g\ 2 represent the atom-atom interactions and are expressed 
as 

2irh 2 a jk 111 

9jk = — , = — + — (39) 

m jk m jk rrij m k 

in terms of the s-wave scattering lengths an and 022 between atoms in the same 
component and a\% between atoms in different components. The s-wave scattering 
length can be tuned by the Feshbach resonance [4~2l H5]. The most striking difference 
between the GP model and the Abelian-Higgs model is that the gauge field is not a 
dynamical variable because the atoms in the condensate are electrically neutral. The 
vector potential A is just an external field and is generated by (i) the rotation of the 
system A = (f2 x r)/2 [52] or (ii) a synthesis of the artificial magnetic field by the 
laser- induced Raman coupling between the internal hyperfine states of the atoms [53] . 

To show the similarity of the formulation between the two-component BECs and 
the field theories in the previous section, we consider the simple case of a homogeneous 
system without a trapping potential Vj = 0. By setting mi = m 2 = m and gn = #22 — g, 
the potential in the hamiltonian Eq. (|38|) can be written as 

-\ 2 



n*) = +G?i2-s)i*ii 2 i^2i 



2 



^ ^W-l^l 2 ), (40) 



2 

where p, = (fi\ + and the constant term has been removed. For g 12 = g and 

A*i = ^2-, the hamiltonian has a global U(2) symmetry. When these conditions are not 
satisfied, the U(2) symmetry is broken to U(l)xU(l). Note that for g 12 < g, the second 
term in Eq. ( 1401) implies that the two components prefer to overlap spatially, while 
for g 12 > g the condensates tend to be segregated to reduce the overlap region. By 
introducing the characteristic length scale 

« = vk' (41) 
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the hamiltonian can be written in a dimensionless form 



H 



d 3 x 



-8li(\^, 



l(V 



1^5 



iA)^\ 2 + - (tfV 
4 v 



V ) 



Ag\V l \ 2 \* 2 \ 



(42) 



with the dimensionless parameters A = 2g//i£ 3 , v 2 = ft£, 3 /g, Ag = (gi 2 — g)/^ 3 fx, and 
5fi = (pi — ix 2 )/jx. For 8/x = 0, Eq.f l42|) is a similar model, discussed by Shifman 
and Yung as a "toy model" The difference from this toy model is the absence of 
the electromagnetic energy, which suggests that the GP model is a counterpart of the 
Abelian-Higgs model in the limit e — > oo but for finite A. 



3.2. Vortices in two- component BECs 

We first briefly describe the vortex states in two-component BECs; see Ref. [20J for the 
details. We confine ourselves to the case Ag = 0, i.e. the SU(2) symmetric case for 
8 fx = and the axisymmetric structure of the vortex states. Since this structure has 
been studied very well, we only consider the analogy with the vortices found in the 
Abelian-two-Higgs model. 

The solution of the vortex state can be obtained from the analysis of the coupled GP 
equations derived from Eq. ( 1421) . The axisymmetric (singly-charged) vortex states are 
characterized by ^(r) = fi(r)e tqid and we consider the most simple case (qi, q 2 ) = (1, 0). 
By applying the boundary condition fi(r = 0) = 0, /i(r = oo) = v, f 2 (r = 0) = 0, 
and f 2 (r = oo) = 0, we can obtain a vortex structure consisting of the circulating ^ 
component that surrounds the non-rotating component at the center, as shown in 
Fig. [71 Since the total density does not vanish at the vortex core, this vortex can be 
called a coreless vortex. If the parameter 5 fx is zero, the system does not allow vortices 
to exist since the vacuum manifold is S 3 which does not have non-contractible loops. 
Thus, the size of the vortex can be arbitrary; the core of the vortex actually fills the 
entire space, in which case the meaning of the vortex is completely lost. Since the gauge 
field is absent in the system, there is no intrinsic mechanism to stabilize the vortex in a 
topological manner like a semilocal vortex in the Abelian-two-Higgs model. Actually, 
an extrinsic mechanism such as the rotation A or the trapping potential Vj allows a 
stable coreless vortex even for Six = [20J . 

When 8 Li > (< 0), the SU(2) symmetry of the system is explicitly broken to U(l) 
x U(l) and the order parameter space is [U(l)x U(l)]/U(l) = U(l) ~ S 1 with the 
ground state = v and j^l = (|^i| = and j^l = v), leading to the formation of 
stable vortex solutions in the order parameter (^2)- However, for smaller values of 
\8ix\, ty 2 (\l/i) can survive only around the vortex core. With increasing \8ix\, the vortex 
core shrinks by an accompanying decrease in the population of ty 2 fifx) a t the core. For 
large values of \8fx\, only the conventional singular-type vortex of remains and the 
field vanishes everywhere in space; the symmetry is restored in the vortex core and 
the core size is fixed by the standard healing length £ = h/ ^/2m/i 1 ( 2 ). 
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Figure 7. Axisymmetric vortex structure of two-component BECs for vortex winding 
number (<?i, (ft) = (1, 0). The parameter values are A = 2 and Ag = in Eq. (|42|). The 
U(2)-symmetry breaking parameter S/j, is chosen to be 0.02, 0.1, and 0.2. (a) Radial 
profile of l^i 1 2 (solid curve) and l^) 2 (dashed curve), (b) Total density p (solid curve) 
and effective magnetic field B z (dashed curve). 



Although the real gauge field is absent, we can define the effective vector potential 
by noting that the gauge field is coupled with the configuration of \I/ in the limit e — > 00 
as seen in Eq. ( TTOj) : 

A cff = -2 ^4 (43) 

In terms of the neutral superfluid, this effective vector potential is equivalent to the 
superfluid velocity field v s . In addition, we can define the effective magnetic field 
through the relation 

B cff = V x A cff = u s , (44) 

which corresponds to the vorticity u s = V x v s . Figure [7] also depicts the distribution 
of B e ff z . In the coreless vortex, the vorticity is continuously distributed around the 
vortex core at which |v| vanishes. As 5fi increases, the vorticity shrinks together with 
the \l/ 2 -component. In the limit of ^2 0, the vortex becomes a singular type and 
the vorticity is concentrated on the origin as a delta function type, whose behavior is 
similar to the B z of the ANO vortex in the e — > 00 limit. 

The behavior of the vorticity distribution can be understood as follows. The velocity 
field of Eq. (j43|) can be written as 

v s = V6 - cos 6Vip, (45) 

which consists of two parts: the current due to the global U(l) phase 6 = 61 + 62 and 
that due to the variation of the polar angle 9 and the azimuthal angle (p — 62 — Q\ of 
the unit sphere (S 2 ). Here, the angles of S" 2 can be related to the pseudospin field 

tfto-tf /^**2 + ^2*i .*** 2 -*2*i l^il 2 - |*2 |2 ~ 



W V l^ll 2 + l^2| 2 ' |^l| 2 + |^ 2 | 2 ' |^l| 2 + |^ 2 | 2 , 

(sin 6 cos (f, sin 9 sin ip, cos 9) . (46) 
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From Eq. (1451) . the vorticity can be written as 



w s = Vx (V6) - cos#V x (V<p) + (V0) x (sin0Vy>) . 



(47) 



The first and second terms in Eq. f|47|) vanish except at the point singularity with 
1^1 = due to the vortex cores. When the vortex core in one component is filled 
by the other component, the singularity of the first and second term is canceled out, 
which allows the non-singular (coreless) vortices described by the third term. The scalar 
product [(V6) x (sin#V0)] • ds is equal to the infinitesimal solid angle covered by the s 
orientations within the infinitesimal plane ds. By using the pseudospin, the third term 
can be rewritten as 



with the Levi-Civita symbol e a ^ 7 . This is known as the Mermin-Ho relation for the s 
texture [54] . 

3.3. D-brane soliton 

Next, we turn to the massive case Ag ^ 0. When the condition Ag > is satisfied, 
the model has two distinct minima as (i) |^i| = v and j^l — 0, (ii) |^i| = and 
|^2| = In the vacuum (i), the field has mass y/ Agv. The fluctuating field 
around the vacuum 5\^i\ = \^\\ — v has mass yXv. In the vacuum (ii), the roles of 
the \I/ fields interchange, as well as their masses. The energies in these two vacua are 
necessarily degenerate because of the Z 2 symmetry <-> apparent in Eq. (j4"2|) . 
which is spontaneously broken. Therefore, there must exist a domain wall interpolating 
between vacua (i) and (ii). 

Let us assume that the wall lies in the z = plane and impose the following 
boundary conditions (\I>i,\I>2) r -> (v,0) T at z -> -00, and (\l>i,\l>2) T ~> (0,w) T at 
z — > 00. The typical profile of the domain wall solution is shown in Fig. [HJ The 
amplitude of ^ decreases as v(l - e^~ 2vz ) as it approaches to the domain wall, while 
that of ^2 increases as e^~^ vz from zero moving away from the domain wall. Thus the 
structure of the domain wall has a two-component structure with thicknesses (y/Xv)^ 1 



and (V Ag v)^ 1 . These two length scales are associated with the total density and 

the pseudospin s z , as seen in the bottom panel of Fig. [8j Generally, the solution of 
the domain wall can be written as \l>i = ^id{z — z^)e~ lLp l 2 and \1> 2 = ^2d{z — z )e tip / 2 , 
containing two parameters — the wall center z and a phase <p = 9 2 —9 1 . Both amplitudes 
of are equal to each other at z = z . The occurrence of z is due to the breaking of 
translational invariance by the given wall solution, while (p is also due to the breaking 
of global U(l); the relative combination of the phase (p is not broken in either vacuum 
(i) or (ii) but is broken along the domain wall, and consequently there appears a U(l) 
Nambu-Goldstone mode localized around the wall. These features satisfy a part of 
the requirement that the domain wall in two-component BECs can be referred to as a 
D-brane soliton |15j . 



-e a ^s a Vs/3 x Vs. 



■7 



(48) 
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Figure 8. Structure of the domain wall in two-component BECs. Top panel shows 
the profile of l^il 2 and l^l 2 , and bottom panel represents the same solution in terms 
of the total density p and the pseudospin s z . 



Since the formulation is similar to the strong coupling limit e — > oo of the Abelian 
massive Higgs model given by Eq. (fl4|) . the corresponding analysis in Sec. 12.31 can 
be directly applied to our problem. According to the Abelian-Higgs model, the vector 
potential and the adjoint field are not independent dynamical valuables, and their forms 
are given by the profile of \1/ through the relations Eq. ( ITOj) and Eq. ( l2~Tj) . For Ag > 
the potential in Eq. ( H2l) can be written as 



V(V) = ^(tfttf - v 2 ) 2 + ^(Eeffl - M) 2 ^ - 5/i(|^i| 2 - |^ 2 | 2 ), (49) 

where M = diag.(M/2, -M/2), M 2 = pAg, and 

M M \^ 1 \ 2 -\^ 2 \ 2 
S c ff = — s z = — . (50) 



2 

The form of Eq. p91) is similar to the potential in Eq. (fl4"|) . but the mass parameter is 
now dependent on the total density M = M(p) [523- Taking this into account, we can 
map the hamiltonian Eq. fj42|) to the generalized massive 0(3) NLcrM |19j : 

~ \ 2 



H = jd 3 x [(V^) 2 + P - $^(Vs a ) 2 + p (v. - A 
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+- A {p-jy + ^p\l-sl)-5ws z _. 



(51) 



Here, we have additional contributions to Eq. ( 122|) . namely, the inhomogeneous effect 
of the total density and the kinetic energy of the superflow velocity (third term in the 
right-hand side of Eq. f|51|) ) . The latter contribution survives in the GP model because 
the vector potential A is just an external field; in the gauge model, A = A = A e ff 
and this term disappears. The density inhomogeneity can be neglected if we consider 
the limit A — > oo, where the radial degree of freedom is frozen as = p = v 2 . In 
this limit the exact solution of the domain wall as well as the D-brane solitons may 
be obtained from Eq. ( I32p . However, we should pay attention to the properties in 
this limit, because the kinetic energy due to the velocity fields v s becomes divergent if 
there are singular vortices in the system. This divergent contribution may be relaxed 
by depleting the total density at the vortex core. Thus, the fixed total density is not a 
good approximation in the BEC system. Detailed analysis shows that, by taking into 
account the density inhomogeneity, the qualitative features of the topological solitons do 
not change, though quantitative profile of the solution is changed through modification 
of the mass parameter M(p) = y/pAg [56J. 

Note that there is a significant difference in the topological structure of the D- 
brane soliton between the GP model and NLaM. In the GP model, the total density p 
decreases exponentially with decreasing z along the vortex core because the nonrotating 
component does not enter energetically into the vortex core [left panel of Fig. [9]^a)] [57] . 
Since the total density vanishes (p — >■ 0) at the singular vortex core, the pesudospin 
is ill-defined there. In this sense, the lump shrinks to a singular vortex for a finite 
distance and thus we can see a boojum connecting a singular vortex to the wall. This 
is in contrast to the case in the NLaM. A lump configuration continues to infinity, 
avoiding the singularity since p is fixed. Hence, a boojum as a singular defect should be 
positioned at z — > — oo in this model [right panel of Fig. [9]^a)]. 

A typical profile of a boojum structure in two-component BECs is shown in the 
bottom of Fig. M, where a single vortex along the z-axis in the ^2 component is attached 
to the domain wall. The numerical solutions show that the total density p remains 
non-zero along the domain wall but vanishes along the vortex core for a finite distance. 
Therefore, we see a boojum connecting a singular vortex to the wall. Above the boojum, 
the vortex becomes a coreless vortex, forming a lump structure in the z =const cross- 
section. 

This statement can be confirmed by considering the distribution of the boojum 
charge density defined in Eq. ( [271) . Since t m of Eq. ( 1271) vanishes for e — > 00, we 
redefine the boojum charge by rescaling as e 2 t m — > t m . Using Eq. (!27j) . we write the 
boojum charge of this system as 



Using Eqs. (fTUl) and (fTTl) in the NLcrM limit to calculate B e fj, we can show that the 
boojum charge density is identically zero. This is because the boojum is relegated to 



V-(S cff B cfr ) = -V-(Ms 2 B eff ). 



(52) 
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Figure 9. Structure of the D-brane soliton in two-component BECs for the simplest 
situation in which a straight vortex is attached to a domain wall, (a) Schematic 
illustration of the D-brane soliton in the NLcrM and the GP model (BEC case). The 
pseudospin on the domain wall has s z = 0. The panels (b), (c), and (d) show the 
numerical solution obtained by solving the GP equation with the parameters Ag = 1 
and Sfi = 0.01718. (b) Profile of the total density, (c) Profile of s z . (d) Vector plot 
of the boojum flux density S ffB c ff in Eq. (|52|) in the small region enclosed by the 
rectangle in (c). The vector only shows the radial and z-component because of the 
cylindrical symmetric structure. The color represents the magnitude of the vector 
Scff-Bcff ■ The mark represents the region where the boojum charge becomes dense. 
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z = ±00. In contrast, the boojum charge in the BEC is localized around the connecting 
point of the coreless vortex and the singular vortex without moving away to z — ±00. 
It is difficult to visualize the boojum charge because it is strongly localized around 
this connecting point, which is needed to calculate the derivative with high accuracy. 
Instead of the boojum charge, in Fig. [9] (d) we show the distribution of the boojum flux 
density S e gB c fj. The flux is concentrated around the connecting point, which implies 
that the boojum is actually localized there as schematically shown in Fig. [9] (d). The 
distribution of the boojum charge in more general configurations of D-brane solitons 
can be studied by full 3D numerical calculation, which will be reported elsewhere. 

4. Conclusion and remarks 

We have shown that an analogue of topological solitons, such as D-brane solitons, 
known in gauge theoretical models can be realized in two-component BECs, because 
the GP model can be regarded as a counterpart of those gauge models. The GP model 
corresponds to the e — > 00 limit of the gauge theory. This fact gives the correspondence 
between two models: (i) The vector potential can be represented as a superfluid velocity 
field, (ii) The adjoint field reduces to the ^-component of the pseudospin density, 
proportional to the density difference l^i] 2 — l^l 2 - We can obtain the analytic solution 
of the D-brane solitons in the strongly interacting limit A — > 00, where the radial 
degree of freedom of ^ is frozen and the model is reduced to the nonlinear a model. 
This approximation can be applied to the BEC problem only in the qualitative level, 
because the density inhomogeneity is crucial around the topological defect. Based on 
the correspondence with the gauge theoretical model, we introduced the boojum charge 
in two-component BECs and showed that the charge is localized around the point where 
a singular vortex terminates on a domain wall. 

We believe that atomic BECs are a promising candidate for demonstrating an 
analog of D-brane physics in a laboratory. The D-brane soliton can be realized as an 
energetically stable solitonic object in phase-separated rotating two-component BECs 
[15] . The long life time of D-brane solitons merits the study of various dynamical 
phenomena, e.g., oscillation modes of strings and branes and nonlinear dynamics 
such as brane-antibrane annihilation, which was proposed as a possible explanation 
for the inflationary universe in string theory. Although brane-antibrane annihilation 
was demonstrated to show the topological defect creation in superfluid 3 He |13j, its 
physical explanation of the creation mechanism of defects is still unclear. In atomic 
BECs, Anderson et al. observed the creation of vortex rings via the dynamical (snake) 
instability of a dark soliton [58J , where the nodal plane of a dark soliton in one component 
was filled with the other component and then the filling component was selectively 
removed with a resonant laser beam. Here, the snake instability can be understood 
as a phenomenon similar to tachyon condensation in quantum field theory [16]. It is 
quite interesting to note the similarity between tachyon condensation in BECs caused 
by domain-wall annihilation and that by brane annihilation in string theory. We also 
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proposed that, when strings are stretched between a brane and an antibrane, namely 
when the filling component has vortices perpendicular to the wall, "cosmic vortons" can 
emerge via a similar instability p2]. All of these phenomena can be monitored directly 
in experiments. 

We do not intend to establish a perfect connection between our D-brane soliton 
in BECs and those in string theory. Strictly speaking, such a perfect connection is 
impossible because the superstring theory is a 10-dimensional theory. Also, our system 
does not contain two important ingredients of original string theory: (i) supersymmetry 
and (ii) relativistic invariance. For (i), the authors in Refs. [3j H] have discussed 
an analogue of D-branes in a field theoretical model, which is similar to our model. 
Although their formulation actually possesses supersymmetry, it does not give an 
essential effect to the solution of the D-brane soliton; even for a supersymmetric theory, 
the solutions are constructed only by the boson part. In this sense, the supersymmetry 
has no effect on the solution. Regarding the correspondence to the original string theory, 
we note that there are interesting proposals [591 160] in which the supersymmetry can be 
included by preparing the Boson-Fermion mixture in an optical lattice. The influences 
of the optical lattice and fermions in our system merit further study. For (ii): the 
relativistic invariance is not significant for our results because we confine ourselves to 
the stationary problem. For the relativistic problem, the time evolution of the scalar 
field is governed by the second-order time derivative, while for the non-relativistic case 
the time derivative is of first order. Therefore, the relativistic invariance is irrelevant 
to our stationary problem. Of course, it is important when we consider the dynamic 
behavior of the D-branes or the strings. However there have been many studies on non- 
relativistic D-branes, applying them to AdS/condensed matter physics correspondence 
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